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Measurement connects the world of quantum phenom- 
ena to the world of classical events. It plays both a pas- 
sive role, observing quantum systems, and an active one, 
preparing quantum states and controlUng them. Surpris- 
ingly - in the light of the central status of measurement 
in quantum mechanics - there is no general recipe for de- 
signing a detector that measures a given observable ^. 
Compounding this, the characterization of existing detec- 
tors is typically based on partial caUbrations or elaborate 
models. Thus, experimental specification (i.e. tomogra- 
phy) of a detector is of fundamental and practical impor- 
tance. Here, we present the reaUzation of quantum de- 
tector tomography [2-4]: we identify the optimal positive- 
operator-valued measure describing the detector, with no 
ancillary assumptions. This result completes the triad, 
state [5-11], process [12-17], and detector tomography, 
required to fully specify an experiment. We character- 
ize an avalanche photodiode and a photon number resolv- 
ing detector capable of detecting up to eight photons 11181 . 
This creates a new set of tools for accurately detecting and 
preparing non-classical fight. 

Von Neumann's postulate of the reduction of the quantum 
state by measurement is now generally accepted to be a lim- 
iting case of a more general theory of quantum measurement. 
However, even within this general theory it is not known how 
to incorporate the complete chain of apparatus components 
in a derivation of the actual measurement: Braginsky wrote, 
"the Schrodinger equation cannot tell us the connection be- 
tween the design of the measuring device and the nature of 
the measurement uJ]." Measurement is increasingly becom- 
ing a driving component in quantum technologies such as 
super-resolution metrology lll9ll . Heisenberg-limited sensitiv- 
ity 1I20I1 . and quantum computing Bill . Input states and dy- 
namical processes are accepted as resources for quantum tech- 
nologies and therefore the techniques of quantum state tomog- 
raphy (QST) [5-11] and quantum process tomography (QPT) 
[12-17] have been developed to measure them. A distinct 
omission is that of the experimental tomography of detectors, 
which would enable more accurate classification of measure- 
ment types, objective comparison of competing devices, and 
precise design of new detectors. This omission is even more 
striking given that the tomography of states and processes are 
predicated on a well characterized detector. In this paper, we 



extend previous theoretical descriptions of detector tomogra- 
phy {l-A} and, by means of efficient numerics based on con- 
vex optimization [22] we characterize two quantum detectors. 
Characterizing a detector consists of determining its corre- 
sponding positive operator valued measure (POVM). Given 
an input state p, the probability p„ p of obtaining detection 
outcome n is 



Pn,p = tr[/9 7r„], 



(1) 



where {7r„} is the detector POVM. In state tomography, an 
unknown p is characterized by performing a set of known 
measurements, each on many identical copies of the state in 
order to estimate j<„. From this estimate one can invert equa- 
tion ^ to find p. The interchangeability of p and 7r„ in equa- 
tion ^ shows that detector tomography plays a dual role to 
state tomography. Now, measuring a set of known probe states 
{p} allows us to characterize an unknown detector, and thus 
find {tt„}. For these operators to describe a physical measure- 
ment apparatus, they must be positive semi-definite, 7r„ > 0, 
and J2„ ^" ~ ^' ensuring positive probabilities that add up 
to one. In addition, the operators {p} must be chosen to be 
tomographically complete, i.e. form a basis for the operator 
space of 7r„. 

In the specific case of optical detectors, lasers provide us 
with an ideal tomographic probe: the coherent state \a). By 
transforming the magnitude \a\ through attenuation (e.g. with 
a beamsplitter) and the phase arg [a] by optical delay, we 
can create a tomographically complete set of probe states 
{|a) (a|} (the existence of the P-function is a proof of com- 
pleteness). Remarkably, with coherent state probes, the mea- 
sured statistics are themselves a full representation of the de- 
tector in the form of the Q-function ||2[], 



Q„{a) = -(a|7r„|a) = -p„,, 

TT TT 



(2) 



Since Q„(a) of each POVM element contains the same in- 
formation as the element 7r„ itself, this is already detector to- 
mography. Predictions of the detection probabilities for ar- 
bitrary input states can then be calculated directly from the 
Q-function representation. Unfortunately, experimental er- 
rors and statistical fluctuations can cause a simple fit to the 
Q-function to be consistent with unphysical POVM elements. 
Due to this we ultimately wish to directly find the POVM el- 



ements {7r„} that are closest to the measured statistics, while 
constraining them to be physical. 

We now turn to the description of the experimental realiza- 
tion, shown in Fig. 2 (see methods). The first detector was 
a commercial single-photon counting module based on a sil- 
icon avalanche photodiode (APD). It has two detection out- 
comes, either outputting an electronic pulse (1 -click) or not 
(0-clicks). Past evaluation of the detector has shown that the 
1 -click outcome is mainly associated with the arrival of one 
or more photons, although dark counts and afterpulsing can 
also create this outcome. The 0-click event is mainly associ- 
ated with vacuum at the input or photons lost due to non-unit 
efficiency of the photodiode. Having only two outcomes, this 
detector cannot directly measure the incoming photon num- 
ber if it is above one. The second detector circumvents this 
by splitting the incoming pulse into many spatially or tempo- 
rally separate bins, making unlikely the presence of more than 
one photon per bin. Subsequently all the bins are detected 
with two APDs. Photon-number resolution results by sum- 
ming the number of 1 -click outcomes from all the bins. This 
time-multiplexed detector (TMD) is not commercially avail- 
able but can be constructed with standard tools lIlSll . Ours has 
eight bins in total (four time bins in each of two output fibres) 
and thus nine outcomes - from zero to eight clicks, making it 
capable of detecting up to eight photons. The added complex- 
ity and greater number of outcomes of this detector provide a 
more challenging test for detector tomography. 

For both detectors we first allowed the phase of a to drift. 
We observed no variation in the outcome frequencies, as ex- 
pected from a detector without a phase-reference. This simpli- 
fies the experimental procedure, requiring us to control only 
the magnitude of a (as has been done for tomography of a 
single photon [23]). A detector with no observed phase de- 
pendence will be described by POVM elements diagonal in 
the number basis. 
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simplifying henceforth the reconstruction of 7r„ . 

For a POVM set {7r„} containing only diagonal matrices 
that are each truncated at a number state M, we can rewrite 
equation (|2]i as a matrix equation. 



p = Fn. 



(4) 



For an N outcome detector, Pu x n contains all the measured 
statistics, F^xm contains the D probe states a, ai^ . . . , an, 
and Haixn contains the unknown POVM set (matrix sub- 
scripts are the matrix dimensions). For a coherent state probe, 
Fi,k — Ickip*^ exp (— |aip)/fc!. This can easily be reformu- 
lated for a probe in a mixed state, as was done to model 
the laser technical noise (see methods). The optimal physi- 
cal POVM consistent with the data can be estimated through 




Coherent state intensity, |ft|' 



FIG. 1: The detector tomography data. The outcome statistics (red 
dots) are measured as a function of the coherent state magnitude \a\ 
and form an estimate of pn, a for each detector outcome n (number 
of clicks). Since they are proportional to the Q-function Qn (a) for 
each outcome, the statistics directly fully characterize the detector. 
The main plot corresponds to time multiplexed detector (TMD) with 
nine outcomes and the inset corresponds to the avalanche photodiode 
(APD). The vertical statistical error is too small to be seen. From the 
reconstructed POVM elements {it„} we generate the corresponding 
probabilities p„, a — {a|7r„|a) (blue curves). 
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FIG. 2: The experimental setup. A half-waveplate (A/2) and Glan- 
Thompson polarizer (?) are used to vary the amplitude of the probe 
coherent state, which is subsequently attenuated by Neutral Density 
Filters (NDF) and coupled into a fibre (FC) (see methods for more 
details). 



the following optimization problem: 



min{||P-Fn||2+5(n)} 



Af-l 



subject to 7r„ > 0, j, ""n = ^' 



(5) 



where the 2-norm of a matrix A is defined as \\A\\2 = 
(Si 7 l^i.j P)"'^^^- Note that we allow for regularization in the 
form of convex quadratic functions g, related to the condi- 
tioning of the problem, which must not depend on the type 
of detector. This is a convex quadratic optimization problem, 
and hence also a semi-definite problem (SDP) |22] which can 
be efficiently solved numerically. Moreover, in this case, there 
exists a dual optimization problem whose solution coincides 
with the original problem. Thus, the dual problem provides a 
certificate of optimality that we use to verify our solution. 

The measured statistics for each detector outcome (i.e. 
number of clicks) are shown in Fig. 1 for the TMD and for the 
APD. The distributions (equivalent to the Q-function (3„(a) 
of the detector) show smooth profiles and distinct photon 
number ranges of sensitivity for increasing number of clicks 
in the detector Fig. 3 shows the diagonals (the off-diagonals 
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FIG. 3: The optimal physical POVMs. We present the diagonals 
of the reconstructed POVMs represented in the photon-number basis 
for (a) the photon-number resolving TMD and (b) the binary APD 
detector. The TMD POVM elements were obtained up to basis state 
|60>(60| (therefore M = 60), but are shown up to |30)(30| for dis- 
play purposes. The APD POVM elements are shown in full. Stacked 
on top of each Qi" we show ^Q^'^'^'^'^' — 6*^ ''° | in yellow, where 
n is the number of clicks, and rec and teo are the reconstructed and 
theoretical diagonals of POVM element, 7r„. The theoretical TMD 
and APD models are described in methods. 



are zero for these phase insensitive detectors) of the POVMs 
that result from optimization of equation ^ (see methods for 
.g(n)). Note that 7r„, being the POVM element for n clicks, 
shows nearly zero amplitude for detecting less than n pho- 
tons, exhibiting essentially no dark counts. Prominent in an 
otherwise smooth distribution, this sharp feature provides the 
detector with its discriminatory power; n chcks guarantees 
there were at least n photons in the input pulse. To assess the 
performance of the tomography we find the difference (yel- 
low bars in Fig. 3) between the estimated POVM elements 
TT™'^ and a previously developed simple theoretical model of a 
TMD, 1:^^° Q (see methods). The fidelity 



^ = tr v/<^<'v^ 



> 98.7% 



for all n, indicating excellent agreement between the two. 

To visualize the action of the detector, in the special case 
of optical detectors one can plot a Wigner function of each of 
the reconstructed POVM elements, Wn (a, a*). The response 
of the detector to an input state with Wigner function W^, is 
proportional to the overlap. 
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We focus on the one cUck Wigner function Wi (a, a*) for the 
APD (Fig. 4a) and the TMD (Fig. 4b). An APD detector 
is sometimes regarded as a 'single photon detector' but here 
we can see the marked difference between the two Wigner 
functions. Instead, it is the TMD that has a fidelity of 98% 
with a single photon (having experienced a 52.2% loss). Con- 
versely, the APD Wigner function extends to a ^ 1 , having 
significant overlap with photon number states ^ 1. Therefore, 
to use an APD as a 'single photon detector' one must make 



FIG. 4: The Wigner functions of the 'one click' detector outcomes. 
From the diagonal elements of tti for the APD (a) and TMD (b) one 
can generate the Wigner function representing their measurement of 
the optical mode. 



the ancillary assumption that the input beam has insignificant 
components containing more than one photon. Despite their 
differences, both Wigner functions have negative values near 
the origin, indicating the absence of a classical optical ana- 
logue. Consequently, these are both fundamentally quantum 
detectors. 

As quantum technologies advance, detectors are becoming 
more complex, making a black-box approach to their charac- 
terization an important tool. Identifying the exact operations 
of detectors will benefit precision tasks, such as state tomog- 
raphy or metrology. By eliminating assumptions, full charac- 
terization enables more flexible design and use of detectors, 
be they noisy, nonlinear, inefficient, or operating outside their 
normal range. With precise characterization we can ask pre- 
cise quantitative questions about our power to prepare non- 
classical states or herald quantum operations |21]. This opens 
a path for the experimental study of yet unexplored concepts 
such as the non-classicality of detectors. For optical detec- 
tor tomography, a promising avenue for research will be to 
transfer well-established techniques from homodyne tomog- 
raphy (e.g. balanced noise-reduction, direct measurement of 
the Wigner function or pattern functions |25]). Now that it is 
well characterized, the photon counter also provides a unique 
tool for performing non-Gaussian operations, which are criti- 
cal for quantum information processing using the electromag- 
netic field as the information carrier 026112711 . As supercon- 
ducting and semiconductor photon number counters are devel- 
oped, tomography could be used as an objective benchmark to 
compare competing devices. Moreover, for one of these pho- 
ton number counters only an incomplete and empirical model 
is available [28], making detector tomography the best option 
to completely determine its action. We expect detector tomog- 
raphy will become the standard for the adequate calibration of 
all measurement and state preparation devices. 
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Experimental setup 

The pulses of a mode-locked laser travel through a half- 
waveplate (A/2) and a Glan-Thompson polarizer (P) with 
which we varied their amplitude a. We subsequently sent the 
pulses through a beamsplitter (BS) (T = 95%). The reflected 
beam travelled through three neutral density (i.e. spectrally 
flat) filters (NDF) before being coupled into a single-mode 
fiber (FC). The attenuation from all elements, the reflection off 
the beamsplitter, each of the filters, and fibre-coupling, were 
measured individually with a cahbrated power meter, result- 
ing in a total attenuation 7. This power meter was then placed 
in the transmission port of the beamsplitter so that the magni- 
tude of a for the probe state in the fibre was found from P, 
the measured time-averaged power and the pulse rate R via 
|a| — ^P\l{2nRhc). For each value of a we recorded the 
number of times each detection outcome occurred in J trials 
(i.e. laser pulses), which provides an estimate of p„ „. 



Source of light and technical noise 

The input states were generated by a mode-locked Ti:Sapph 
laser with center wavelength A and a FWHM bandwidth 
of AA specifically chosen for each detector. It was cavity 
dumped to reduce its repetition rate R in order to be compati- 
ble with tested detectors. Long term drift of the intensity over 
1 million pulses was < 0.5%. To characterize it, a NIST cal- 
ibrated Coherent FieldMaxII-TO power meter was used (sys- 
tematic error of 5%). In the case of the APD detector (a Perkin 
Elmer SPCM-AQR-13-FC) we set A = 780± 1 nm, AA = 20 
nm, and chose the appropriate rate R = 1.4975 ± 0.0005 
kHz, J = 1472967, and 7 = (5.66 ± 0.08) x 10"^. For 
the TMD detector we set A = 789 ± 1 nm, AA = 26 
nm, R = 76.169 ± 0.001 kHz, J = 38084, and 7 = 
(8.51 ± 0.11) X 10^^. We now evaluate the importance to 
our tomography of the technical noise found at some level in 
all lasers. Our laser randomly varies in energy between sub- 
sequent pulses with a standard deviation of 1.88% ± 0.02% 
of I a I . Attenuated to the signal photon level, as in this ex- 
periment, one might expect the inherent large fractional un- 
certainty in the coherent state to render this technical noise 
insignificant. We test this expectation by modelling the pulse 
distribution as a Gaussian /^(/S) = e"('^"")''/(2'^')/(crV27r) 
centered around a in phase space, with a variance approxi- 
mately equal to that measured, a^ = 0.0004 \a\ . Each probe 
state is then best described by a mixture of coherent states. 
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The detection probability for outcome n is then 
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Comparing our analysis done with pure input states |a) (a| to 
that done with mixed states p^^,^ we find the difference be- 
tween the POVMs obtained was negligible. For example 
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and the largest relative difference between any two 6*^. com- 
ing from a mixed state or a pure state derivation was 1.3%. 
Furthermore the reconstructed probability distributions are so 
close that they are indistinguishable on the scale of Fig. 1. 
This reinforces our earlier expectation that technical noise in 
the laser will be negligible when using single-photon-level co- 
herent states. This differs from homodyne tomography where 
technical noise can shift a strong local oscillator to a nearly 
orthogonal state. 



Discussion of regularization 

Care has to be taken that the optimization problem is well 
conditioned in order to find the true POVM of the detector In 
finding the number basis representation we are deconvolving 
a coherent state from our statistics, which is intrinsically an 
ill-conditioned problem. Similar issues of conditioning have 
been discussed in the context of state and process tomogra- 
phy, see e.g. Refs. 11291 BOll . Due to a large ratio between 
the largest and smallest singular values of the matrices defin- 
ing the quadratic problem, small fluctuations in the probabil- 
ity distribution can result in large variations for the recon- 
structed POVM. This can result in operators that closely ap- 
proximate the outcome statistics and yet contain errant spikes 
in their distribution in photon-number To suppress this ef- 
fect, we penalize the difference ^J." — 0^^-^ (independent of 
the shape of the POVM) by using the regularization g — yS 

with S = Efe,„[^i"^ - ^i+i]^- This is motivated by the fact 
that any realistic detector will have a finite efficiency 77, which 



l,m=0 



necessitates a smooth 9j^ distribution: if G{r) is the proba- 
bility of registering r photons and H{q) is the probability that 
q were present then. 

Consequently, if 6*^: 7^ then fl^+i, 9k+2 etc. cannot be zero, 
but will follow some smooth distribution. Since we do not 
assume any knowledge about the precise loss of our detector 



we simply choose an arbitrary value for y. Varying y by three 
orders of magnitude hardly affects the exact value of the es- 
timated POVM, changing it by only 10%. Furthermore, the 
regularization g ^ y S also proves to be robust to noise up to 
S = 0.2 (varying a -^ a(l + S) across {|Q;)(a|} with a Gaus- 
sian distribution for S). This shows that the regularization's 
main effect is to suppress the ill-conditioning and noise while 
leaving the POVM fitting unaffected. 



The theoretical model of the detector 

The detector tomography does not make use of any phys- 
ical model on the functioning of the detector To verify the 
success of this approach we have compared the outcome of 
the estimation with those POVM elements obtained from a 
theoretical model of the APD and TMD |,24J. The APD is 
treated as a binary detector with a loss of 43.2%. The theo- 
retical TMD assumes: No dark counts, three sequential beam 
splitters with experimentally inferred reflectivities, 50.18%, 
50.60%, and 41.92%, and an overall loss of 52.2% (that best 
fits the data), followed by two perfect APDs. 
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